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Abstract 

We present exact solutions for the zero-temperature partition function of the g-state 
Potts antiferromagnet (equivalently, the chromatic polynomial P) on tube sections of the 
simple cubic lattice of fixed transverse size x Ly and arbitrarily great length Lz, for sizes 
3 and 2x4 and boundary conditions (a) {FBCx, FBCy, FBCz) and (b) 
{PBCx, FBCy, FBCz), where FBC (PBC) denote free (periodic) boundary conditions. In 
the limit of infinite-length, — oo, we calculate the resultant ground state degeneracy 
per site W (= exponent of the ground-state entropy). Generalizing q from Z+ to C, we 
determine the analytic structure of W and the related singular locus B which is the continuous 
accumulation set of zeros of the chromatic polynomial. For the — > oo limit of a given 
family of lattice sections, W is analytic for real q down to a value qc- We determine the 
values of qc for the lattice sections considered and address the question of the value of qc for 
a (i-dimensional Cartesian lattice. Analogous results are presented for a tube of arbitrarily 
great length whose transverse cross section is formed from the complete bipartite graph 
K 
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1 Introduction 



The g-state Potts antiferromagnet (PAF) exhibits nonzero ground state entropy, 5*0 > 
(without frustration) for sufficiently large g on a given lattice A or, more generally, on a 
graph G. This is equivalent to a ground state degeneracy per site W > 1, since 5*0 = 
kB In W . Such nonzero ground state entropy is important as an exception to the third law of 
thermodynamics [0, ^. One physical example is provided by ice, for which the residual molar 
entropy is 5*0 = 0.82 ± 0.05 cal/(K-mole), i.e., Sq/R = 0.41 ± 0.03, where R = NAvogkn §• 
Indeed, residual entropy at low temperatures has been observed in a number of molecular 
crystals, including nitrous oxide, NO and FCIO3 (a comprehensive review is given in Ref. 0). 
In these physical examples, the entropy occurs without frustration, i.e., the configurational 
energy can be minimized, just as in the Potts antiferromagnet for sufficiently large q. 

There is a close connection with graph theory here, since the zero-temperature parti- 
tion function of the above-mentioned g-state Potts antiferromagnet on a graph G = (V, E) 
satisfies 

ZiG,q,T = 0)pAF = PiG,q) (1.1) 

where G is defined by its set of vertices V and edges E and P{G, q) is the chromatic polyno- 
mial expressing the number of ways of coloring the vertices of G with g colors such that no 
two adjacent vertices have the same color (for reviews, see [0-|llO|). The minimum number 
of colors necessary for such a coloring of G is called the chromatic number, x{G). Thus 

W{{G},q) = lim P(G,g)i/" (1.2) 

where n = |V^| is the number of vertices of G and we denote the formal infinite-length 
limit of strip graphs of type G as {G} = lim„_+oo G- At certain special points g^ (typically 
g^ = 0, 1, .., x{G)), one has the noncommutativity of limits |]11 



lim lim P(G,g)^/" ^ lim lim P(G,g)^/" (1.3) 

q—>qs n— >oo n— >oo q^qs 

and hence it is necessary to specify the order of the limits in the definition of W{{G},qs). 
Denoting Wgn and Wnq as the functions defined by the different order of limits on the left and 
right-hand sides of (|1.3|), we take W = Wgn here; this has the advantage of removing certain 
isolated discontinuities that are present in Wnq. Using the expression for P{G,q), one can 
generalize g from Z_|_ to C The zeros of P{G, q) in the complex g plane are called chromatic 
zeros; a subset of these may form an accumulation set in the n ^ 00 limit, denoted i3, 
which is the continuous locus of points where W{{G},q) is nonanalytic. For some families 
of graphs B may be null, and W may also be nonanalytic at certain discrete points. The 
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maximal region in the complex q plane to which one can analytically continue the function 
W{{G}, q) from physical values where there is nonzero ground state entropy is denoted R^. 
The ground state degeneracy per site is an analytic function of real q from large 

values down to the value qc, which is the maximal value where B intersects the (positive) real 
axis. For some families of graphs, B does not cross or intersect the real q axis; in these cases, 
no qc is defined. However, even in cases where no such intersection occurs, B often includes 
complex-conjugate arcs with endpoints close to the positive real axis, and hence, in these 
cases, it can be useful to define a quantity {qc)eS equal to the real part of the endpoints. We 
shall use this definition here. 

In this work we present exact solutions for chromatic polynomials P{G, q) for sections 
of the simple cubic (sc) lattice with fixed transverse size x Ly and arbitrarily great 
length Lz, for cross sections x Ly ~ 3 x 2 and 4x2. These calculations are carried out 
for the cases (a) {FBC^, FBCy, FBC^) (rectangular sohd) and (b) {PBC^, FBCy, FBC,) 
(homeomorphic to an annular cylindrical solid), where FBCi and PBCi denote free and 
periodic boundary conditions in the i'th direction, respectively. We shall use the notation 
{Li)-p and (ivi)p to denote free and periodic boundary conditions in the i'th direction, so 
that, for example, the 3 x 2 x L;^ sections of the simple cubic lattice with the boundary 
conditions of type (a) and (b) are denoted Sp x 2f x {Lz)-p and 3p x 2f x {Lz)-p, respectively. 
For each family of graphs, taking the infinite-length limit — > oo, we calculate ^^({(7}, g), 
B, and hence qc- 

We also present corresponding results for a tube of arbitrarily great length whose trans- 
verse cross section is formed from the complete bipartite graph Km,m, for the cases m = 2 
and 3. Here the complete graph Kn is defined as the graph consisting of n vertices such that 
each vertex is connected by edges (bonds) to every other vertex, and the complete bipartite 
graph Km,n is defined as the join + K^, where the join of two graphs G and H, denoted 
G + H is the graph obtained by joining each of the vertices of G to each of the vertices of H. 
Although the complete bipartite graphs are not regular lattices as studied in statistical me- 
chanics, they are useful since they allow us to obtain exact results for cases of high effective 
coordination number. 

There are several motivations for this work. We have mentioned the basic importance 
of nonzero ground state entropy in statistical mechanics and physical examples of this phe- 
nomenon. From the point of view of rigorous statistical mechanics, exact analytic solutions 
are always valuable since they complement results from approximate series expansions and 
numerical methods. We have defined the point qc above in terms of the function 1^({G}, q). 
This point has another important physical significance: for the n — > oo limit of a given family 
of graphs, {G}, the g-state Potts antiferromagnet has no finite-temperature phase transition 
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but is disordered for all T > if g > qc{{G}), and has a zero-temperature critical point for 
Q = Qc{{G}) [0. For the Potts model on the (infinite) square lattice, via a mapping to a 
vertex model, it has been concluded that q = 3 [|15|. However, the value of Qc is not known 
for any lattice of dimension three or higher. One of the main motivations for our study is the 
insight that one gains concerning the dependence of Qc on the coordination number A = 2d 
of a dimensional Cartesian lattice E*^ for the case d = 3. Furthermore, although infinite- 
length sections of higher- dimensional lattices with fixed finite {d — l)-dimensional volume 
transverse to the direction in which the length goes to infinity are quasi-one-dimensional 
systems and hence (for finite-range spin-spin interactions) do not have finite-temperature 
phase transitions, their zero-temperature critical points are of interest. Finally, in addition 
to the physics motivations, the present results are of interest in mathematical graph theory. 



Some related earlier work is in Refs. |T3|1- [M - 

Our exact calculations of Qc for the infinite-length limits of tube sections of the simple 
cubic lattice yield information relevant to estimates of Qc for the infinite cubic lattice. The 
point here is that as the area of the transverse cross section of the tube increases to infinity, 
the corresponding sequence of exact qc or (gc)efr values is expected to converge to a limit. 
For the tube sections considered here, which have free longitudinal boundary conditions, this 
limit is a lower bound for the true Qc of the infinite simple cubic lattice. The reason that one 
can only say that it is a lower bound is that for these families with free longitudinal boundary 
conditions the respective limiting curve B exhibits a complex-conjugate pair of prongs that 
protrude to the right. It is possible that, as the area of the transverse cross section goes to 
infinity, the endpoints of these prongs will extend over and meet on the real axis, thereby 
defining a point qc which could lie to the right of the limit of the Qc points for each of the tubes 
with finite transverse cross section. In this case, limL^^Ly^oo Qcisc, {Lx)bcx x {Ly)BCy x oop) 
is not equal to, but instead less than, the value qc{sc) for the infinite simple cubic lattice. 
Indeed, for finite-width, infinite-length strips of the triangular lattice with free longitudinal 
{z) boundary conditions and periodic transverse y boundary conditions, ImiLy^oo Qc{{Ly)p x 
oop) ^ 3.81967... while gc = 4 for the infinite triangular lattice [^. This occurs because 
of the above-mentioned phenomenon in which the endpoints of the complex-conjugate pair 
of prongs protruding to the right on the limiting curve B move inward toward the real axis 
and, as — > oo, finally touch each other, thereby defining the true gc = 4 about 5 % greater 
than the above limit at ~ 3.82. The presence of these types of prongs protruding to the 
right on the limiting curves was also observed for a variety of infinite-length, finite-width 2D 
lattice strips with free longitudinal boundary conditions (for both free and periodic transverse 
boundary conditions) |]3D|, ^ |5D|. Calculations for Lp x oop (i.e. cylindrical) 



strips of the square lattice have been carried out in [^, ^ with widths L extending 
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up to 13 |5^, with qc{sq, 13p x oop) — 2.916. This value is within about 3 % of the value for 
the infinite square lattice, qdsq) = 3. As the width increases, the endpoints of the prongs 
do move in toward the real axis. These calculations for Lp x oop strips of the square lattice 
are consistent with either of the two possibilities, that linii^oo qd^q, Lp x oop) is equal to, 
or slightly less than, qdsq). 

It is interesting that the exact calculations of the singular loci B for infinite-length, finite- 
width strips of various lattices with periodic longitudinal boundary conditions |^ (and either 
free or periodic transverse boundary conditions) yielded loci without such complex-conjugate 
prongs (or line segments on the real axis) [Tl|, |36|, ^ ^ . For these families of strips 
it was found that B always crossed the positive real axis at g = and at a maximal point, so 
qc is always defined; furthermore, Lp x oop, the value of qc is a monotonically nondecreasing 
function of the width. For Lp x oop strips, although q^ is not a nondecreasing function of 
the width, it is, for a given width, closer to the value for the infinite 2D lattice if one uses 
periodic, rather than free transverse boundary conditions, as one would expect since the 
former minimize finite-size effects. 

Some definitions from graph theory will be useful here. The degree of a vertex in a 
graph G = (V, E) is the number of edges to which it is attached. A graph G defined to be 
A-regular iff each of its vertices has the same degree, A. This is, in particular, true of an 
infinite regular lattice, where this degree is the coordination number. Even if some vertices 
have different degrees from others, one can define an effective coordination number, 

AefT = lim 4^ . (1.4) 

IvHoo \V\ 

For example, for a strip of the square lattice of size of the type {L;j.)f x {Ly)F, in the limit 
in which the length Lx c>o with fixed finite width, Ly, 

2 

Aefr(sg, {Lx)f x (Ly)p; ^ oo) = 4 - — . (1.5) 

Ly 

For the same strip, but with periodic transverse boundary conditions, if Lj^ > 3 to avoid 
double edges, we have 

Acfr(sg, (L^)p x {Ly)p ; ^ oo) = 4 . (1.6) 

For the {Lx)p x {Ly)p x {Lz)p section of the simple cubic lattice with free boundary conditions 
in all three directions, the interior vertices have degree 6, the vertices on the surface away 
from the corners have degree 4, and the corner- vertices have degree 3. For the minimal case 
2p X 2p X {Lz)-p., there are no interior vertices in the 3D sense, so that 

Aeff(sc, 2f X 2f X (L,)f; ^ cx)) = 4 . (1.7) 
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For the next case {Lx)f x 2^? x {Lz)-p with > 3 there are also no interior vertices in the 
3D sense; in this case 

Assise, (L,.)f X 2f X (L,)f; FBC,, FBCy, FBC,; ^ oo) = 5 - ^ . (1.8) 

For the ^ oo hmit of the x Ly x section with L^; > 3, L^, > 3, we have 



Aeff(sc, (L,)f X (L,)f X (L,)f; ^ oo) = 6 - 2^^|±^ . (1.9) 

We shall also use the corresponding formulas when periodic rather than free boundary condi- 
tions are imposed in one of the transverse directions. Of course, in the usual thermodynamic 
limit of the Cartesian lattice E'^, with Lj — > oo, i = 1, .., d and lim^.^oo Li/Lj equal to a finite 
nonzero constant, the effective degree is Aeff(E'^) = 2d independent of boundary conditions. 
In this case, although free boundary conditions in one or more directions entails vertices 
of lower degree than 2d, these vertices occupy a vanishing fraction of the total vertices in 
the thermodynamic limit. Concerning the chromatic number, for the sections of 3D lattices 
considered here, we have 

X{sc, (L,)f X {Ly)F X (L,)f) = 2 (1.10) 

i.e., these graphs are bipartite. For the sections of the simple cubic lattice with periodic 
boundary conditions in one of the transverse direction, taking this to be the x direction with 
no loss of generality, 

, . ^ X . ^ X / r \ \ f 2 if Lj,. is even /i i . \ 

X(.c, (L.)P X (L,)f X (L.)f) = I 3 if^^ i^^^d (1.11) 

For an infinite regular lattice having coordination number A, a rigorous upper bound on 
Qc was derived in using the Dobrushin theorem: 

ge < 2A . (1.12) 

Aside from the general property that, for cases where there is a Qc, this point satisfies Qc > 2, 
we are not aware of a published lower bound on Qc- 

We have noted that in the infinite-length limit of given family of lattice strip graphs, B 
does not necessarily intersect or cross the real q axis, and for such a family, there is no qc- It 
has been found that for a graph consisting of a strip of a regular lattice, a sufficient condition 
for B to cross the (positive) real q axis is that there be periodic boundary conditions in the 
longitudinal direction, i.e., the direction in which the length goes to infinity as ^ oo 
lTT|-i4l. 
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An interesting question is how qc depends on the lattice coordination number, equivalent 
to the vertex degree for an infinite lattice. For families of x Ly lattice strip graphs with 
periodic longitudinal {Lx) boundary conditions, whose ^ oo limits are guaranteed to 
have a Qc-, it is found that the value of Qc is a nondecreasing function of the effective vertex 
degree Agfr- Thus, for the family of cyclic strips of the square lattice with free transverse 
boundary conditions, besides the the li? x oo^ case with A = 2 and qc = 2, one finds the 
following results: (i) 2i7 x oo^? (A = 3) yields gc = 2 |Tl|, (ii) Sp y< oop (Aeg = 10/3) yields 
~ 2.3365 H, (in) 3^ x oop {A^s = 7/2) yields ~ 2.4928 [||]. When one makes both 



the longitudinal and transverse boundary conditions periodic, the vertex degree is fixed. In 
this case, Qc does not necessarily increase with increasing width, Ly. For example, for strips 
of the square lattice with toroidal boundary conditions and hence A = 4, the 3p x oop strip 



yields gc = 3 while the 4p x oop strips yields the smaller value Qc — 2.78 P9 



When one considers strips without periodic longitudinal boundary conditions, there may 
or may not be a qc- We have found that for strips with free transverse and longitudinal 
boundary conditions, qc or (gc)efT, where the latter can be defined, is a monotonically in- 
creasing function of Ly, and since in these families of lattice strips, A^s is a monotonically 
increasing function of Ly, this means that {qc)cs is also a monotonically increasing function 
of Aeff. Thus, for Ip x oop with Agg = 2 and 2p x oop with A^s = 3, S = and there 
is no (gc)cfr5 while 3p x oop with Acg = 10/3 yields qc = '2 |3^, 4p x cxop with A^s = 7/2 



yields ^ 2.3014 g, 5p x cx)p with A^s = 3.60 yields g^ ^ 2.428 [||, etc. up to 



8p X oop with Acff = 3.75, which yields arc endpoints at g ~ 2.6603 ± 0.0013z and hence 
(gjeff = 2.6603 0. 

In contrast, for strips with free longitudinal and periodic transverse boundary conditions, 
Aeff is fixed, e.g., equal to 4 for strips of the square lattice. Insofar as one expects qc or (gc)eff 
to depend on Agg, it is not clear, a priori, how these quantities should behave as functions of 
Ly. Indeed, one finds that qc or (gc)efr have no monotonic dependence on Ly. For example. 



for 3p X oop, i3 = 0, 4p X oop yields g^ ~ 2.3517 0, 5p x cx)p yields g^ ~ 2.6917 pl 
but 6p X oop yields qc — 2.6132 For 7p x oop, i3 has arc endpoints near the real axis 

at g ^ 2.7515 ± 0.0025^, so (gc)cff ^ 2.7515 g^]. For the cases {Ly)p x oo up to L = 11, it 
is found empirically that qc or (gc)efr is monotonically increasing separately for the even-Lj^ 
and odd-Ly sequences |^ . 

Since qc has thermodynamic significance as the value of g above which the Potts antifer- 
romagnet has no finite-temperature phase transition and is disordered even at T = 0, the 
value of qc for a lattice with dimensionality d > 2 should be independent of the boundary 
conditions used to take the thermodynamic limit. (Strictly speaking, this is not true for 
d = 1, since no qc is defined for the infinite-length limit of the line graph, while qc = 2 for the 
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infinite-length limit of the circuit graph, but this may be considered to be an exception due 
to the special nature of this graph.) For the infinite 2D lattices where Qc is known exactly, it 
is also a monotonically increasing function of the lattice coordination number. Specifically, 



as noted above, on the square lattice (A = 4), one has = 3 [16|. For the kagome lattice. 



again with A = 4, one finds = 3 |T31, and on the triangular lattice with A = 6, one has 



qc = ^m 



These exact results are all consistent with the inference that the value of qc increases as 
a function of the coordination number for the thermodynamic limit of a regular lattice or, 
more generally, for infinite-length strips or tubes of regular lattices with prescribed boundary 
conditions in the transverse directions. A plausibility argument for this is as follows. For 
q > qc, the zero-temperature Potts antiferromagnet has nonzero ground state entropy per 
site, S = ksliaW > 0, i.e., ground state degeneracy per site W > 1. This entropy reflects the 
fact that the number of ways of carrying out a proper coloring of the lattice or lattice strip 
increases exponentially with the number of vertices. Roughly speaking, the constraint that 
no two adjacent vertices can be assigned the same color is more restrictive the greater the 
number of neighboring vertices there are, i.e., the greater A^s is. Therefore, as one increases 
AefT for a fixed q, W decreases This is borne out by exact and numerical calculations of W (for 
2D lattices, see, e.g.. Fig. 5 of [|lT| or Fig. 6 of |2^). In the context of d-dimensional Cartesian 



lattices, this means that q^ should increase as a function of d. Such a monotonicity relation 
is valuable to have since the value of qc is not known exactly for lattices with dimension 
d> 3. 



It is useful to define a ratio of the actual value of qc to the upper bound in ( 1.12 ), namely, 



= ^ . (1.13) 

Evidently, Rq^ = 0.5 for the infinite-length limit of the circuit graph and Rq^ = 0.375 for 
the infinite square lattice. Similarly, for other exactly known cases, the actual value of qc 
is considerably less than the rigorous upper bound (|1.12|) . This suggests that it may be 
possible to improve the upper bound (|1.12|) . 

A generic form for chromatic polynomials for a strip graph of type Gs, or, more generally, 
a recursive family of graphs composed of m repetitions of a basic subgraph, is |2^ 



P{Gs,m,q) = J2 CG.,m^GsM)r (1-14) 

where cc^jiq) and the Ng^,\ terms (eigenvalues) Xcsjil) depend on the type of strip graph 
Gs including the boundary conditions but are independent of m. Here, Nx < dim(T) , where 
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dim(T) is the dimension of the transfer matrix in the Fortuin-Kasteleyn representation, and 
the difference, A^oa = dim(T) — Nx, is the number of zero amphtudes (coefficients). The 
coefficients cg^j can be regarded as the multiphcities of the eigenvalues XgsJ (Jot some 
real positive q values, these coefficients can be zero or negative, so that this interpretation 
presumes a sufficiently large real positive q, followed by analytic continuation to other values 
of q). 

2 Tubes of the Simple Cubic Lattice with Free Trans- 
verse Boundary Conditions 

The computations of the transfer matrix and the chromatic polynomials were performed 
following the methods described in Section 3 of Ref. The idea is to construct the 



partition function by building the lattice layer by layer. In all the computations we have 
chosen the Fortuin-Kasteleyn representation of the transfer matrix. Thus, our basis will be 
the connectivities of the top layer, whose basis elements v-p are indexed by partitions V of 
the single-layer vertex set V^. We shall abbreviate delta functions S{ai, a^) as ^1^3. Moreover, 
we shall also abbreviate partitions V by writing instead the corresponding product of 
delta functions: e.g., in place of V = {{1,3}, {2,4}, {5}} we shall write simply V = 61^362^4,- 
The transfer matrix can be expressed as 

T = VH (2.1) 
where the matrices H and V are defined as 

H = n l^-J--'] (2-2) 
V = n M + ^J • (2-3) 



xgyo 

In these formulae is the single-layer edge set, and J^^x' and are respectively the "join" 
and "detach" operators whose action on the elements of the basis is as follows 

Jxx'^V = ^V»xx' (2.4) 

Dx^v = I ^ I (2.5) 

' \ gv-p if{2;} eV ^ ' 

where V • xx' is the partition obtained from V by amalgamating the blocks containing x and 
x' (if they were not already in the same block); and V\x is the partition obtained from V 
by detaching x from its block (and thus making it what we term a "singleton"). In eq. (|2.3| ) 
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we have written the formulae for H and V in the zero-temperature Potts antiferromagnet 



hmit; the general expressions can be found in [47 



Finally, the partition function can be written as 

Z{L, xLyX L,- g) = ■ H (VH)^^-i ■ Vid (2.6) 

where "id" denotes the partition in which each site x G is a singleton, and u""" is defined 
by 

. VP = gl^l . (2.7) 

In the zero-temperature limit the horizontal operator H is a projection. This implies that 
we can work in its image subspace by using the modified transfer matrix T' = HVH in place 
of T = VH, and using the basis vectors 

= H vp (2.8) 

in place of v-p. Note that w-p = if P contains nay pair of nearest neighbors in the same 
block. In the following subsections we will list the basis P = {vp}, although we performed 
the actual computations with the basis {wp}. 

The computation of the limiting curves B was performed with either the resultant and 
the direct-search methods (See ref. for details). The computation of the zeros of the 
chromatic polynomials was done using the package MPSolve designed by Bini and Fiorentino 



5^ , |53| , whose performance is superior to that of Mathematica for this particular task [^ |. 



2.1 2f X 2f X {Lz)y Section of the Simple Cubic Lattice 

Before proceeding to our new results, we note the identity 

sc(2f X 2f X (L,)f) = sg(4p x (L,)f) . (2.1.1) 

The left-hand side of this identity is evidently a tube of the simple cubic lattice with a 
minimal-size transverse cross section, viz., a single square; the right-hand side of the identity 
is a strip of the square lattice of width 4 squares, and periodic transverse boundary conditions. 
The chromatic polynomials for the family of graphs on the right-hand side of this identity 
were calculated in section 7 of [pT| and, in the infinite-length limit, W and B were determined 



(see Fig. 3(a) of [0). In Ref. [0, the longitudinal direction was taken as L^, while we 
take it to be here. In ||3y]-|3^, a generating function formalism was developed and used; 
for a given strip graph Gg of length m, the chromatic polynomial P{{Gs)m, q) is given as the 
coefficient of z^, where z is the auxiliary expansion variable: 

oo 

r(G„g,^) = 5^P((G,)„,g)z™, (2.1.2) 

m=0 
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This generating function is a rational function of q and z. For the strip considered here, of 
length Lz = m edges, 

V{Gs,q,z) 



ao + aiz 



where 



1 + hiz + 62-2^ 



«o = q{q - - 3g + 3) 



(2.1.3) 



ai 



= -q{q - 1)(/ - + 16g^ - 13g + 5) 
hi = -g^ + 8q^ - 29q'^ + 55g - 46 



(2.1.4) 
(2.1.5) 
(2.1.6) 
(2.1.7) 

(Here we have converted the numerator coefficients from to be in accord with the labelling 
convention of eq. ( |2.1.2| ); in the labelling convention of |^T|, the strip with length m + 1 
edges.) The chromatic polynomial can equivalently be written as in |2.6| : in the basis P = 
{1, (5i,4 + 62,3} it is given by 



q^ - Vli + 61g^ 



169g^ + 269g^ - 231g + 85 



Z(2f X 2f X (L,)F;g) 



g(g-l)(g2-3g + 3) 
2q{q-\f 



T(2f X 2f 



(2.1.8) 



The transfer matrix T(2f x 2f) is equal to 

41 



r(2F X 2f) 



51g + 28g2 - 
-5 + 2g 



g,3 _^ ^4 



-12 + 14g - 6g2 + q^) 
5 - 4g + g2 



(2.1.9) 



Note that in general, we should have considered an additional element in the basis, namely 
the partition 51^453^2; however, its amplitude vanishes identically. This can be easily under- 
stood by noting that this particular graph ( |2.1.1[ ) is planar and 5i,453^2 is a crossing partition. 

For this family, B consists of two complex-conjugate arcs, a self-conjugate arc that crosses 
the real axis at g ^ 2.3026, and a line segment on this axis extending from q ~ 2.2534 to 
q ^ 2.3517, which latter point is thus qc. It should be noted also that although this graph is 
included here as the minimal-transverse-size member of a family of x Ly x sections of 
the simple cubic lattice, it is degenerate in the sense that it is actually planar, in contrast to 
the others that we shall study here. There are eight endpoints in B] two of them are given by 
the real values q ~ 2.3026 and q ~ 2.3517; the other six endpoints are the complex-conjugate 
pairs q ~ 0.7098 ± 2.0427 i, q ~ 1.9923 ± 1.5942 and q ~ 2.9953 ± 1.4266 i 



2.2 3f X 2f X (L^)f Section of Simple Cubic Lattice 

The section of the simple cubic lattice with the next larger transverse cross section is the 
3f X 2f X {L2)yi or equivalently, 2f x 3f x [Lz)yi family. For this family, the dimension of 
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the transfer matrix is 13. (In general the transfer matrix has dimension 15, but two basis 
elements, namely 51,552,6^3,4 + 52,4^3, 5^1,6 and 5i,3,552,4,6) can be eliminated as they have zero 
amplitudes.) Since the entries in the transfer matrix are rather long, we relegate them to 
the Mathematica file transf er_sc .m that is available with this paper in the LASL cond-mat 
archive. We have computed the chromatic zeros for = 15 and 30, i.e., n = 90 and 180 
vertices, respectively. These are shown in Fig. ^ This value of is sufficiently large that 
these zeros give a reasonably accurate indication of the location of the asymptotic limiting 
curves comprising the locus B. As is generally true, for the larger value of L^, and hence n, 
the chromatic zeros move slightly outward, approaching the limiting curve B from smaller 
values of \q\. 



Zeros sc lattice = 3p Ly = 2p 




Re(q) 

Figure 1: Chromatic zeros for the 3f x 2f x {Lz)f section of the simple cubic lattice, for (a) Lz = 15, i.e., 
71 = 90 (□), (b) Lz = 30, i.e., n = 180 (o). 

The limiting curve B consists of two complex-conjugate arcs, a self-conjugate arc that 
crosses the real g-axis at g ^ 2.8645, and a horizontal real segment from q ^ 2.8566 to 
q ^ 2.8723. This latter value corresponds to the value of qc for this family. This information 
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Table 1: Results for tube graphs of the simple cubic lattice. The table shows the relation between the 
degree (coordination number) A of a A-regular family or the effective degree Aes, and qc, if a qc exists, for 
the infinite-length limit of the family. Boundary conditions are indicated with a subscript F for free, P for 
periodic, in a given direction. In column q^ an asterisk indicates that B does not actually cross the real 
axis, so that, strictly speaking, no qc is defined, but arcs on B end very close to the real axis, at the positions 
given. Similar results are listed for the tubes with Km.m cross sections. 



Gs 


A 


Aeff 


Qc 


ref. 


sc, 2f X 2f X oop 




4 


2.3517 








sc, 3f X 2f X oof 




4.33 


2.8723 


here 


sc, 4f X 2f X oof 




4.50 


3.1498 ±0.0021* 


here 


sc, 2p X 2p X oof 




4 


2.3517 





sc, 3p X 2p X oop 




5 


3.3255 ± 0.0184 


here 


sc, 4p X 2p X oop 




5 


3.3623 ±0.006H* 


here 


K2,2 X OOp 




4 


2.3517 








-f^3,3 X OOp 




5 


3.0452 ± 0.0082 r 


here 


^^4^4 X CXDp 




6 


3.6743 ± 0.0085 r 


here 



is hsted, together with that from other lattices, in Table |1|. There are 16 endpoints: two are 
the real values q ~ 2.8566 and q ~ 2.8723; the other 14 are the complex-conjugated pairs: 
q ~ 0.2943 ± 1.9150 i, q ~ 1.9371 ± 2.4040 i, q ^ 2.6374 ± 1.0937i, q ~ 2.8250 ± 1.4087^, 
q ^ 2.6255 ± 1.4647 i, q ~ 3.2410 ± 1.2920 i, and q ~ 2.2645 ± 2.4605 i Finally, there are 
four T points at g ~ 2.670 ± 1.117^ and q ~ 2.689 ± 1.500 z. 

2.3 4f X 2f X (L^)f Section of Simple Cubic Lattice 

The section of the simple cubic lattice with the next larger transverse cross section is the 
4f X 2f X {Lz)-F-, or equivalently, 2f x 4f x {Lz)-f-, family. For this family, the dimension 
of the transfer matrix is 156. However, there are 20 basis elements that correspond to 
zero amplitudes, so the effective dimension of the transfer matrix is 136. This transfer 
matrix T(4f x 2f), as well as the vectors v and Ujd can be found in the Mathematica file 
transf er_sc .m. We have computed the chromatic zeros for = 20 and 40, i.e., n = 160 
and 320 vertices, respectively. These are shown in Fig. ^ 

The limiting curve B consists of eight complex-conjugate arcs; none of them crosses the 
real axis, so, strictly speaking, there is no value of qc- However, we can define (gc)eff ~ 3.1498. 
This value is larger than the corresponding value for the family 3f x 2f x Lp qc ~ 2.8723 
and larger than the value for the square lattice qdsq) = 3. There are 11 pairs of complex- 
conjugate endpoints: q ~ 0.073 ± 1.747 i, q ~ 1.582 ± 2.845 i, q ~ 1.684 ± 2.871 i, q ~ 



12 



Zeros sc lattice = 4p Ly = 2p 




Re(q) 



Figure 2: Chromatic zeros for the 4f x 2f x {Lz)f sectfon of the simple cubic lattice, for (a) — 20, i.e., 
n = 160 (□), (b) = 40, i.e., n = 320 (o). 
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2.379±2.503i, q ~ 2.513±1.569 q ~ 2.661 ±1.816 i, q - 2.697±1.8892, q ~ 3.034±1.3672, 
q ~ 3.035±1.385 i,q — 3.150±0.002 i, and q ~ 3.383±1.157i. Finally, there are three pairs of 
complex-conjugate T points: q ~ 2.381 ±2.404i, q ~ 2.695± 1.793i, and q ~ 2.972 ± 1.578 i 

3 Tubes of the Simple Cubic Lattice with Transverse 
Periodic Boundary Conditions 

The method we have used to compute the chromatic polynomials and the transfer matrix is 
the same as in Section ^. The only difference is that we enlarge the single-layer edge set . 

3.1 2p X 2p X {Lz)-p Section of Simple Cubic Lattice 

This family has (trivially) the same chromatic polynomials as the family sc{2-p x 2f x {Lz)y- 
(This is not true for the finite-temperature Potts model partition function, or equivalently, 
the Tutte polynomial; however we only deal with the chromatic polynomial here.) Thus, the 
transfer matrix and chromatic polynomials are given respectively by (|2.1| ) and (|2.6|) . 

3.2 3p X 2p X {Lz)-p Section of Simple Cubic Lattice 

In general, finite-size artifacts are minimized if one uses periodic boundary conditions in as 
many directions as possible. Hence, it is useful to calculate P and B for the same section 
of the simple cubic lattice as in the previous section, but with periodic boundary conditions 
imposed on the longer of the two transverse directions. This family is again bipartite. In 
this case, the dimension of the transfer matrix is 4. 

In the basis P = {1, ^2,4 + ^2,6 + ^1,5 + (^3,5 + + ^3,4, ^1,5^2,6 + ^2,4^3,5 + ^3,4^1,5 + ^3,4^2,6 + 

^i,6'^3,5 + ^1, 6(^2,4, ^1,6^3,4 + ^1, 5^2,4 + 52fi5z,b} ^ the transfer matrix can be written as 




T(3p X 2p) 



1 2(-3 + g) 12-6g + g2 _2(-4 + g) 



(3.2.1) 



V 



-2 -4(-3 + q) 15 - 7g + g2 y 



where 



Tl4 



Til 



13 



12 



1089 - 1578g + 1054g^ - 417g^ + 103g^ - 15g^ + q^ 
6(-292 + 394g - 231^^ + 7Aq^ - 13g^ + q^) 
6(110 - 117g + 51g2 - HgS + q^) 
3(154 - 149g + mq^ - I2q^ + q^) 



(3.2.2) 
(3.2.3) 
(3.2.4) 
(3.2.5) 
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T21 



22 



-90 + 71g - 20g^ + 2q^ 
203 - 189g + 7lq^ - 13g^ + 



(3.2.6) 
(3.2.7) 
(3.2.8) 
(3.2.9) 



-109 + SOg - 21g2 + 2q^ 
-81 + 52g - 12g2 + g3 




2 + g)(-l + g)g(-13 + lAq - 6g2 + g3) \ 
6(-2 + g)3(-l +g)g 
6(-2 + g)2(-l +g)g 
3(-3 + g)(-2 + g)(-l + g)g / 



V 



(3.2.10) 



v 



1 








(3.2.11) 



v 



We remark that in this case the most general basis contains an additional partition: (5i,552,6'^3,4+ 
(^2,4^3, 5(^1, 6- This one can be dropped, as it corresponds to a vanishing amplitude. 

We have computed the chromatic zeros for = 15 and 30, i.e., n = 90 and 180, 
respectively. These are shown in Fig. |^. The limiting curve B contains three pairs of self- 
conjugated arcs. Although the arcs do not cross the real axis, so that, strictly speaking, no 
qc is defined, (gc)eff — 3.33, which is larger than the value gc = 3 for the square lattice. The 
locus B has 12 endpoints: q ^ 0.5061 ±2.6413 q ~ 2.1301 ± 2.4407i, q ~ 3.0251 ± 2.5249 z, 
q ~ 3.0412 ± 1.2643 z, q ^ 2.9328 ± 1.1238 z, and q ~ 3.3255 ± 0.01839 z. 

3.3 4p X 2p X {Lz)-p Section of Simple Cubic Lattice 

In this case we have 56 basis elements. However, there are 11 trivial basis element that lead 
to a vanishing amplitude. The transfer matrix and the vectors v and Ujd are listed in the 
Mathematica file transf er_sc .m. Among the other 45 basis elements, we have numerical 
indications that there are 19 additional vanishing amplitudes. 

The limiting curve B (see Fig. ^ contains eight self-conjugate arcs. As in the 4p x 2p x oop 
case, this locus does not actually cross the positive real q axis, but again has endpoints that 
lie very close to this axis, and we obtain (gc)eff — 3.36. As expected, this is larger than 
the value (gc)eff — 3.33 that we found for the tube 3p x 2p x oop. There are nine pairs of 
complex-conjugate endpoints: q ^ —0.095 ± 2.700 z, g ^ 1.803 ± 3.031 g ^ 2.176 ± 3.196 z, 
g ~ 2.841 ± 1.954i, g ~ 2.910 ± 2.080i, g - 3.263 ± 1.333i, g - 3.279 ± 1.377i, g ~ 
3.362 ± 0.006 z, and g ^ 3.892 ± 1.542 z. Finally, there is one pair of complex-conjugate T 
points: g ~ 3.204 ± 1.665 i. 
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Zeros sc lattice = 3p Ly = 2p 




Re(q) 



Figure 3: Chromatic zeros for the 2p x 3p x {Lz)^ sectfon of the shiiple cubic lattice, for (a) — 15, i.e., 
n = 90 (□), (b) = 30, i.e., n = 180 (o). 
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Another property that was observed in our earher calculations of chromatic polynomials 
for strips of 2D lattices is that, for a given type of transverse boundary conditions, for 
infinite-length strips with free longitudinal boundary conditions, as the width increases, the 
number of arcs on B increases and the endpoints of these arcs move in such a manner as 
to reduce the gaps between the arcs. Here we observe the same qualitative behavior as the 
area of the transverse cross section of the tube section increases. If this trend continues for 
progressively larger transverse widths or cross sectional areas, then the number of arcs could 
increase without bound as one approaches the respective infinite 2D or 3D lattices. Now in 
the two cases where the loci B have been calculated exactly for regular lattices, namely the 
d = 1 lattice (with periodic boundary conditions; for free boundary conditions, B = 0) and 
the triangular lattice (defined as the limit L oo of Lp x oop strips this locus has no 
prongs or endpoints. Thus, one could imagine that as the width or cross sectional area of the 
strips or tubes increases to infinity, the endpoints of arcs join so that the gaps between these 
arcs disappear, and prongs either have lengths that go to zero or have endpoints that join to 
form closed boundaries. We also observe that, for a given width or, for the present families 
of tube sections, for a given transverse cross sectional area, there tend to be fewer arcs when 
one uses periodic rather than free boundary conditions for the transverse direction(s). This 
is in accord with the fact that calculations for strips with periodic longitudinal boundary 



conditions [0, 0, ^ ^ ^ found no prongs (or line segments) on the respective loci 
B, i.e. in all cases studied, these loci did not contain endpoints. Finally, our calculations 
are consistent with the expectation that as the area of the transverse cross section goes to 
infinity, the outer envelope of the locus B approaches a limit, which crosses the real axis at 
q = 0, qc{sc), and other point(s) between these two. 

4 Further Remarks on the Structure of B 

In the introductory section, we discussed some differences in the locus B that depend on 
whether one uses periodic or free longitudinal boundary conditions. Here we include some 
further remarks relevant to our present results. It has been established that the structure of 
B inside of its outer envelope and, in particular, the question of whether and where it crosses 
the positive real axis between g = and q = qc{{G}) are dependent upon the boundary 
conditions used. For example, exact calculations for Lp x oop and Lp x oop strips of the 



triangular lattice show that the respective loci B cross the real axis at g = 2 [^, ^ 
corresponding to the fact that the Ising antiferromagnet has a (frustrated) T = critical 
point on these strips; however the loci B obtained in |23, 31 1 for Lp x oop strips or in [3^ for 



Lp X oop strips of the triangular lattice do not, in general, cross the real axis at g = 2, nor is 
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Zeros sc lattice = 4p Ly = 2] 




Re(q) 



Figure 4: Chromatic zeros for the 2p x 4p x {Lz)^ section of the simple cubic lattice, for (a) — 20, i.e., 
n = 160 (□), (b) = 40, i.e., n = 320 (o). 
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this crossing obtained for the infinite-width hmit of the cyhndrical strips calculated in . 
Similarly, for Lp x cxop and Lp x oop strips of the square lattice it was found that B passes 
through q = 0, q = 2, and a maximal value, q = qc{{G}) [|ll|, ^ In contrast, for the 
corresponding strips of the square lattice with free longitudinal boundary conditions, it has 
been found that B does not contain g = or, in general, q = 2, and while the arc endpoints 
nearest to the origin move toward this point as the width increases (leading to the inference 
that for infinite width, B would pass through q = 0), there is no analogous tendency of arcs 
on B to elongate toward q = 2 0, |3l|, ^ |5^. There is also no tendency for the arcs on 
the loci B for our tube sections of the simple cubic lattice with free longitudinal boundary 
conditions to move toward g = 2 as the cross sectional area increases. 

A second point concerns the positions of the leftmost arcs. We find from the exact 
calculations reported here that for the infinite-length tube sections of the simple cubic lattice 
of the form {Lx)p x {Ly)p X(yop with sufficiently large transverse cross section, namely L^, = 4, 
Ly = 2, B contains support in the Re{q) < half-plane. This suggests that for this family of 
tube sections of the simple cubic lattice with periodic transverse boundary conditions, as the 
transverse area L^Ly oo, the complex-conjugate curves on B could approach the origin 
from the Re{q) < half-plane, as happens for the limit L ^ oo of Lp x oop strips of the 
triangular lattice and sufficiently wide strips of the square lattice of the form Lp x oop 



43], Lp X oop PI, 13, 13, and Lp x oop [||, ||, ||] 



5 Behavior of Qc for with large d 



For the (i-dimensional Cartesian lattice E , Mattis suggested the ansatz |24 



iy(E",g)^l + i^^. (5.1) 

This agrees with the general results 1^({C}, q) = q — 1 for d = 1, and, for the square lattice 
with g = 3 yields the estimate W{sq,q = 3) = 3/2, which is within 4 % of the known 
result W{E^,q = 3) = (4/3)^/2 = 1.5396.. pf. Mattis addressed the question of 4(g), 
i.e. the lower critical dimensionality of the g-state Potts antiferromagnet, below which it is 
disordered for T > 0. An equation yielding 4 as a function of g can also be solved to yield 
gc as a function of d, so dc{q) and qdd) constitute equivalent information about the system. 
Mattis argued that dc{q) could be estimated by noting that W is a measure of disorder, and 
if it is significantly greater than 1, then the system would be sufficiently disordered that one 
would not expect there to be a phase transition at T > 0. Taking the criterion that W < 2 
as the demarcation value for which a zero-temperature phase transition could occur, this 
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yields the result 

= i^i^ . (5.2) 

As noted, this may equivalently be regarded as an equation for as a function of d and for 
d = 3, this ansatz gives 

gc(rf = 3)~4.15 (5.3) 

Our exact results W and qc for tube sections of the simple cubic are consistent with this 
estimate from the ansatz ( ^.l] ). 

It is also of interest to ask what the behavior of qc is for large d on Cartesian lattices. 
The ansatz (|5.1| ) yields the asymptotic behavior 

qc ~ rr~l f*^^ d oo . (5.4) 
met 

This is in agreement with the upper bound ( |1.12| ) 

Qc < 4c? (5.5) 
and evidently is a smaller and smaller fraction of this upper bound as d gets large, with 



6 Family {K^^ 



m, 



It is also useful to calculate B and study the dependence of qc on vertex degree for infinite- 
length limits of other families of tube graphs. We report calculations here for a family of 
tube graphs whose transverse cross section is the complete bipartite graph Km,m.- The graph 
K„i,m is A-regular graphs with A = m. We construct our tubes with the Km,m transverse 
cross section connected lengthwise Lz times, so that each vertex of one Km,m subgraph is 
connected "vertically" to the corresponding vertex of the next Km,m- This recursive family 
of graphs is denoted as (-ft'^.m)'^"'- The value of Aeg is for this family 

AefT {{K^,n,f-,FBC./, ^ oo) = m + 2 . (6.1) 

The computational method is the same as in the simple-cubic families: the transfer matrix 



and the partition function are computed using the same formulae (|2.1| ) / (|2^ ). The only 
difference is the single- layer edge set E^. 
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6.1 Family {K2,2] 



The family is trivially equivalent to the family sq(4p x Lp), so by it is 

also equivalent to the family sc(2f x 2f x {Lz)f) The transfer matrix, written in the basis 
P = {1, 61^2 + ^3,4}^ takes the same form as (|2.1.9| ); and the partition function is equivalent 
to (|X8|). 



6.2 Family (Ks^s)^^ 

The transfer matrix for the family (^^3 3)'^^ has dimension 5. In the basis P = {1, ^1,2 + ^1,3 + 

^2,3 + (^4,5 + (^4,6 + (^5,6, (^1,2,3 + (^4,5,6, ^1,2^4,5 + 51,2(^4,6 + (^1,2(55,6 + ^1^364^5 + (5i,354,6 + 5l,3^5,6 + 
^2,3^4,5 + ^2,3^4,6 + (^2,3^5,6, '^1,2,3'^4,5 + '5l,2,3'^4,6 + '^1,2,3'^5,6 + ^,5,6^1,2 + (^4,5,6^1,3 + ^,5,6^2,3}, We 

can write the transfer matrix as 



where 



TiK: 



33) 





T12 


^13 


Tl4 


Tl5 


\ 


T21 


^22 


^23 


^24 


T25 




T31 


T32 


^33 


^34 






1 


2(-3 + g) 





T44 


2(-2 + g) 




V 


1 





-3 


2-q 


J 



(6.2.1) 



Til 


= 1234 - 1747g + 1137^^ - 437^^ - 


h 105g^ - 15q^ + q^ 


(6.2.2 


T12 


= 6(-252 + 337g - IQSg^ + 65g^ - 


I2q^ + g') 


(6.2.3 


Tl3 


= 2(36 - 56g + 33g2-9g^ + g^) 




(6.2.4 


Tl4 


= 9(89 - 94g + 43g2-10g3 + g^) 




(6.2.5 


Tl5 


= 6(-25 + 24g - + g^) 




(6.2.6 


T21 


= -113 + 91g-27g^ + 3g^ 




(6.2.7 


T22 


= 174 - 153g + 54g2 - lOg^ + g4 




(6.2.8 


T23 


= -10 + 13g - 6q^ + 




(6.2.9 


T24 


= 3(-40 + 28g - + g^) 




(6.2.10 


T25 


= 29-21g + 4g2 




(6.2.11 


T31 


= 55 - 27g + 3g2 




(6.2.12 


T32 


= -3(29- 17g + 3g2) 




(6.2.13 


T33 


= 9-12g + 6g^-g^ 




(6.2.14 


T34 


= -18(-3 + g) 




(6.2.15 


T35 


= -3(7-5g + g2) 




(6.2.16 



21 



= 12-5g + g^ 



(6.2.17) 



Finally, 



/ (-1 + g)g(31 - 47g + 28g2 - 8q^ + q^) \ 



\ 



6(-l + g)g(-7 + lOg - Sg^ + 
2{-l +qfq 
9(-l + g)g(3-3g + g2) 

6(-l + g)2g 



(6.2.18) 



/ 1 \ 








(6.2.19) 



We remark that in this case there is an additional element of the basis 5i,2,3'^4,5,6 that 
should be taken into account in general. However, it corresponds to a vanishing amplitude. 

Chromatic zeros for (fsTs.s)^^ with = 15 and = 30 are shown in Fig. ^, as well as the 
limiting curve B. The limiting curve B contains six connected pieces. None of them crosses 
the real axis. Thus, strictly speaking, there is no qc defined. However, by extrapolating the 
closest points to the real axis we get qc — 3.045, which is slightly greater than the value for 
the square lattice qc = 3. 

There are 14 endpoints: q ~ 0.8197±2.9764 z, q ~ 1.9761±2.5559z, q ~ 2.8190±1.5587z, 
q ~ 2.9364±2.5742i, q ~ 3.6220±2.0051 i, q ~ 3.0283±1.3476 i, and g ^ 3.0452±0.008246z. 
There are two complex-conjugated T points at g ^ 2.949 ± 1.870 z. 

6.3 Family {K^,^)^^ 

In this case the transfer matrix has 15 elements. However, three of them correspond to null 
amplitudes, so we have an effective 12- dimensional transfer matrix. This matrix is listed in 
the Mathematica file transf er_Knn_tube .m that is available with this paper in the LASL 
cond-mat archive. 

There are eight connected pieces (See Figure P), and none of them crosses the real axis. 
The closest points to that axis are the complex-conjugated pair g ^ 3.6743 ± 0.0085 z. There 
are ten endpoints (that were computed using the resultant method): g ^ 1.0084 ± 3.7740 
g ^ 1.9104 ± 3.4341 i, q ^ 2.9457 ± 3.2436 i, q ^ 3.0385 ± 2.8658 i, q ^ 3.5456 ± 1.3512 i, 
q ^ 3.6006 ± 3.2332 i, q ^ 3.6260 ± 1.4516 i, q ^ 3.6743 ± 0.0085 i, q ^ 3.7857 ± 2.3839 i, and 
g ^ 3.8460 ± 2.7980 i. There are four T-points at g ^ 3.070 ± 2.904 i, and g ^ 3.567 ± 3.158 i 
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Zeros (K lattice n= 3 




Re(q) 



Figure 5: Chromatic zeros for the (i^a.a)™ graph for (a) m = = 15, i.e., n = 90 (□), (b) m ^ = 30, 
i.e. n = 180 (o). 
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Zeros (K lattice n= 4 




Re(q) 



Figure 6: Chromatic zeros for the {K4,4)"^ graph for (a) ^ 40, i.e., n = 320 (□), (b) m ^ ^ 80, 

i.e. n = 640 (o). 
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7 Conclusions 



In this paper wc have reported exact solutions for the zero-temperature partition function of 
the g-state Potts antiferromagnet on tubes of the simple cubic lattice with various transverse 
cross sections and boundary conditions and with arbitrarily great length. We have used these 
to calculate, in the infinite-length limit, the resultant ground state degeneracy per site W 
and the singular locus B which is the continuous accumulation set of the chromatic zeros. 
In particular, we have calculated the value of Qc or {qc)eS for these infinite-length tubes. Our 
results show quantitatively how this quantity increases as the effective coordination number 
for a given family of graphs increases and are a step toward determining Qc is for the infinite 
simple cubic lattice. We have also presented similar calculations for another interesting 
family of tube graphs whose transverse cross section is formed from the complete bipartite 
graph Km ,m- 
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H = — J^^jj^ ^o-i,cTj with J < 0, which obviously has a Gibbs measure. 
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